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For  the  four  years  of  Stanford’s  research  activity 
under  this  grant,  it  is  believed  that  the  most  significant 
"products"  are  the  substantial  number  of  publications  and 
the  advanced  graduate  students  whose  theses  formed  the  foun¬ 
dations  of  those  publications.  Nearly  all  of  these  students, 
along  with  a  host  of  others  whose  work  was  supported  by  OSR 
over  a  continuous  period  beginning  in  1953,  are  now  con¬ 
structively  employed  in  academia,  industry  or  the  laboratoi — 
ies  of  allied  governments  around  the  world.  Their  names  are 
recorded  as  authors  or  coauthors  of  archival  papers,  SUDAAR 
reports,  reports  of  the  MIT  Aeroelastic  &  Structure  Research 
Laboratory,  etc. 

At  the  beginning  of  the  attached  list  of  references, 
the  Principal  Investigator  has  attempted,  in  roughly  chron¬ 
ological  order,  to  summarize  most  of  the  papers  whose  con¬ 
tents  were  wholly  or  partially  supported  by  the  grant.  Some 
of  these  have  been  published,  in  whole  or  part,  by  archive 
journals  subsequent  to  issuance  of  the  cited  report?  others 
will  be  in  the  near  future. 

Many  opportunities  have  occurred,  and  will  continue  to 
occur,  for  less  formal  communication  of  recent  research  dis¬ 
coveries.  Several  of  these  have  already  been  described  to 
OSR,  for  example  in  the  rejected  proposal  Aero  No.  1-89  sub¬ 
mitted  in  Sept.  1988  and  in  the  annual  Interim  Scientific 
Reports.  As  part  of  the  process  of  completing  his  doctoral 
requirements,  candidate  M.  Ameen  Jarrah  summarized  his  ex¬ 
perimental  program  at  a  Stanford  Fluid  Mechanics  Seminar  in 
December  1988.  By  invitation,  the  Principal  Investigator 
gave  talks  on  the  unsteady  flow,  agile-aircraft  maneuvers 
and  loads  findings  to  engineers  of  Boeing  Commercial  Air¬ 
planes  on  March  88,  1989,  to  Boeing  Military  Airplanes  in 
Wichita  on  June  88,  1989,  and  to  a  seminar  audience  at  San 
Diego  State  University  on  March  8,  1989. 
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As  a  final  observation,  it  is  noted  that  aerospace 
organizations  around  the  world  have  not  all  lost  sight  of 
the  Principal  Investigator’s  career  of  contributions  to 
research  and  teaching  in  unsteady  aerodynamics,  aeroelas- 
ticity  and  related  fields.  In  October,  1987,  he  was  awarded 
the  Ludwig  Prandtl  Ring  by  DGLR,  the  West  Germany  aerospace 
professional  society  (one  of  a  total  of  five  such  recogni¬ 
tions  for  living  Americans).  In  December,  1988,  the  British 
Royal  Aeronautical  Society  selected  him  as  one  of  its  two 
Honorary  Fellows  for  that  year.  In  June,  1989,  the  West  Gei — 
man  research  institution  DLR  invited  him  to  give  the  keynote 
aftei — dinner  speech  on  applied  optimization  for  their  Semi¬ 
nar  on  Optimization  in  Bonn;  this  was  one  of  their  annual 
series  of  high-level  technical  seminars  on  chosen  topics  in 
the  field. 


Summary  of  Research  Prior  to  Mid-1988 

This  activity  has  been  fully  described  in  three  of  the 
aforementioned  Interim  Scientific  Reports,  dated  in  April  or 
early  May  of  1986,  1987  and  1988.  By  way  of  summary  and 
prior  to  the  work  of  Dr.  Jarrah  described  in  more  detail 
below,  it  is  believed  that  the  principal  contributions  sup¬ 
ported  by  the  grant  are  those  extensively  reported  in  the 
dissertations  of  Dr.  van  Niekerk  (published  in  Refs.  I  and 
II),  Dr.  Brandao  (Ref.  VI)  and  Dr.  Azevedo  (Ref.  VII).  It 
should  be  mentioned  that  the  last  two  individuals  were  Bra¬ 
zilian  nationals  and  that,  in  considerable  part,  their  work 
and  attendance  at  Stanford  were  funded  by  that  government. 
Especially  in  the  case  of  Dr.  Azevedo,  however,  there  was 
substantial  involvement  of  the  grant;  this  is  being  recog¬ 
nized  in  the  resulting  publications  in  the  usual  way. 

Azevedo’s  accomplishments  are  regarded  as  particular¬ 
ly  outstanding  and  of  substantial  interest  to  U.  S.  Air 
Force,  with  the  impending  payload  launches  by  the  Titan  IV 
series  of  booster  configurations.  Starting  from  first  prin¬ 
ciples  but  some  25  years  after  the  incidents  which  he  ana¬ 
lyzed,  he  was  able  to  predict  successfully  an  instability  of 
"hammerhead"  payloads  on  ballistic  launch  vehicles.  In  so 
doing  he  coupled  a  linear-elastic  representation  of  the  LV, 
based  on  superposition  of  it  first  three  natural  modes  of 
free-free  bending  vibration,  with  a  transonic,  unsteady  CFD 
code  employing  approximate  Navier-Stokes  equations  with  a 
modified  Baldwin-Lomax  relation  between  shear  stresses  and 
rates  of  strain.  The  latter  was  adapted  from  axisymmetr ic , 
steady-flow  codes  developed  by  Pulliam.  For  an  Atlas-Able 
vehicle  which  encountered  difficulty  of  this  kind  in  the 
early  1960’s,  he  predicted  a  high— q,  transonic  instability 
of  what  aeroelasticians  call  the  "single-degree-of -freedom" 
variety  for  the  17-Hz  second  mode.  The  first  and  third  modes 
were  found  to  be  quite  stable,  a  prediction  which  agreed 
with  in-flight  observation  as  well  as  could  be  ascertained. 
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Summary  of  Recent  Research  in  1988  and  1959 

The  bulk  of  this  report  deals  with  progress*  prior  to 
and  during  the  cited  period*  on  Dr.  Jarrah’s  unsteady  high- 
angle-of-attack  testing  and  preliminary  attempts  to  make 
applications  of  his  results.  The  experimental  effort  relied 
completely  in  the  availability  of  time  in  one  of  the  7'  x 
10-’  low— speed  wind  tunnels  at  NASA  Ames  Research  Center.  In 
this  connection  the  project  was  indebted  since  early  1987  to 
Mr.  Richard  Margason*  Chief*  Fixed-Wing  Aerodynamics  Branch* 
as  well  as  to  Mr.  Tim  Naumowicz*  who  is  an  engineer  assig¬ 
ned  to  that  branch.  They  arranged  for  tests  to  take  place 
during  two  extended  periods*  the  first  in  late  August  and 
September*  1987,  and  the  second  in  April  and  May  of  1988. 
NASA  also  furnished  the  strain-gauge  balance*  laser  illumi¬ 
nation*  large  quantities  of  hardware  and  software  for  data 
handling*  and  personnel  support  before  and  during  the  tunnel 
entries.  The  dollar  value  of  this  support  is  estimated  at 
well  over  $10*000. 

In  November*  1987,  Mr.  Margason  suggested  that  parti¬ 
cular  measurements  of  interest  to  NASA  might  be  included  in 
the  program  through  a  joint  research  interchange  with  Stan¬ 
ford  under  what  is  called  the  NASA  Ames  University  Consoi — 
tium.  The  conversations  resulted  in  the  award  of  Contract 
NCA2-S87,  entitled  "Unsteady  Flow  Measurements  on  Delta  Wing 
Models  Forced  to  High  Angles  of  Attack."  With  Margason  and 
the  Principal  Investigator  as  collaborators*  this  was  for  a 
period  of  one  year  commencing  January  1,  1988,  and  funded  at 
$25,000.  The  collaboration  was  entirely  complementary  to  the 
activity  supported  by  OSR*  and  the  Project  Monitor  was  noti¬ 
fied  in  a  timely  fashion.  More  detail  is  given  about  this 
arrangement  in  the  1988  Interim  Scientific  Report.  The  in¬ 
tention*  if  Dr.  Jarrah  is  able  to  complete  plans  to  visit 
Stanford  during  the  summer  of  1989,  is  that  a  final  report 
on  the  NASA  contract  will  be  issued  in  the  form  of  a  Tech¬ 
nical  Note  presenting  all  of  the  data  obtained  under  the 
program. 

The  OSR-suppor ted  inves-  *tion  during  this  recent 
period  is  believed  to  be  so  we i *  described  in  the  last  pa¬ 
per  prepared  by  Dr.  Jarrah  and  the  Principal  Investigator 
<Ref.  X)  that  an  Appendix  is  attached  hereto  adapted  from 
that  document.  Any  question  about  those  results  or  other 
accomplishments  under  the  grant  can  be  directed  to  the  Prin¬ 
cipal  Investigator,  telephone  (415)  723  4136. 
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IMPACT  OF  FLOW  UNSTEADINESS  ON  MANEUVERS  AMD 
LOAOS  OF  AGILE  AIRCRAFT 
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and 

Holt  Ashley** 

Department  of  Aeronautics  anc  Astronaut  3 cs 
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Abstract 

The  paper  begins  by  reviewing  a  new 
p.-ojram  of  unsteady  airload  measurements 
executed  on  a  family  of  low-aspect-rat id 
delta  wings  and  motivated  by  recent  in¬ 
terest  in  " super maneuver *"  a*  a  capabil- 
: tv  of  the  next  generation  of  aircraft 
designed  for  ai’ — superiority  missions. 
In  transient  pitch  motions  for  time  con- 
staita  and  maximum  o's  which  reproduce 
full-scale,  normal  force  and  other  loads 
signi f icaitly  exceed  steady— state  veluea 
when  x  is  increasing  but  fall  far  below 
on  the  downstioke.  For  a  series  of  exam¬ 
ples  involving  "generic"  super maneuvers 
taken  from  the  literature,  implications 
of  this  discovery  are  illustrated.  Turn 
rates  are  achieved  which  can  be  consi-i 
derably  greater  than  what  one  would  pre-- 
diet  with  steady  wind-tunnel  data.  This 
increase  in  agility  does  not,  however, 
necessarily  require  a  penalty  in  terms 
of  inermased  structural  loada.  A  sim¬ 
plified  "theory*  is  proposed,  trying  to 
show  how  tha  important  influence  of  the, 
leading-edge  vortex  instability  eight 
empirically  be  incorporated  into  load-, 
ing  estimates  for  wings  with  sharpened 
edges.  Conclusions  are  stated  regarding 
the  introduction  of  these  findings  into  . 
thr  design  of  such  aircraft. 


Nomenclature 

AR  <-b“/S>  Wing  aspect  ratio 

b  Wing  span 

c  Wing  chord 

“  Mean  aerodynamic  chord 

Co,  Ct.  Cm.  Cm,  Cw  Coefficients  of 

drag,  lift,  pitching  moment ,  normal 
force  and  rolling  moment  (standard 
aeronautical  definitions) 

Co«  Drag  coeff  icleflT  of  plate  normal  to 
airf low 

Cmw  Pitching  moment  about  model  axis 
D  Drag  force 

g  Acceleration  of  gravity 

•■Research  Assistant)  Ph.C.  awarded  Janu¬ 
ary  1989.  Student  member,  AI AA 

**Pnfessor.  Honorary  Fellow.  AIAA. 


h  Altitude 

k  (-ftc./av)  Reduced  frequency 
K  <*an«ac./SV)  Pitch-rate  parameter 

L  Lift  force 

m  Maas  of  flight  vehicle 

Mw  Pitching  moment  about  pitch  axis 
n*  Normal  load  factor 

N  Normal  force  on  wing 

q  <“trV“/2)  Flight  dynamic  pressure 

0  Angular  velocity  in  pitch 
Re  Reynolds  no.  based  on  midspan  chord 
S  Plan  area  of  wing 

t  time 

T  Thrust  force 

V  speed  of  flow  or  flight 

x,y,z  Ca-tesi an  coordinates  (x  measured 
aft  from  vertex  of  delta  wing) 
x„(t)  x-coordinate  for  vortex  breakdown 
x„  x-coordinate  of  pitch  axis 
X  Azimuth  angle  of  flight  path 


oi  Angle  of  attack 

‘ 6  Angle  of  sideslip 

r  Flight-path  .ingle  above  horizontal 

fl  Sweep  angle  of  wing  leading  edge 
o  density  of  air 

p  Sank  angle  about  velocity  vector 

C  Circular  frequency  of  sinusoid 

Sufrfqr 

<*)  Time  derivative 

<  >.  Cartesian  coordinates  in  horizontal 
plane 

(  )m«  Maximum  value  of  time  function 

(  >»  Reference  value  !c.  is  midspan 
chord  of  wing) 


IfliLttdMsliga 

The  tactical  advantages  of  "super¬ 
maneuvers"  in  short-range  combat  between 
air-superiority  aircraft  were  first 
pointed  out  in  the  open  literature  by 
Herbst  iRefs.  1,  S  and  several  more 
recent  publ ications) .  They  are 
effective  at  very  low  speeds,  where 
transients  of  angle  of  attack  a  to  30“ 
and  above  can  ba  performed  without 
exceeding  acceleration-  tolerances  of 
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aircrew  or  structure.  They  have  also 
motivated  a  groat  deal  of  research  in 
vehicle  dynaaicsi  trajectory 

optimization,  aerodynamic*  and  related 
fields.  Important  studies  of  optimized 
supermarteuver*  by  Well*  Faber  and  Ber¬ 
ger4**  appeared  in  the  early  1900'*.  The 
many  collection*  of  paper*  on  the  topic 
are  typified  by  AGARD  Conference*  <e.g., 
Dietz  and  Due**’)  and  the  USAF  Technical 
Specialists  Meeting4**. 

It  i*  noteworthy  that  published  an¬ 
alyse*  of  high-a  tactic*  have,  for  lack 
of  better  information,  been  forced  to 
rely  on  steady-state  airload  data.  Good 
current  examples  Mould  be  presentations 
(e . g . ,  Anderson4**)  and  a  panel  discus¬ 
sion  at  the  1989  Aerospace  Science*  Mee¬ 
ting.  Similarly,  the  minimum-time  turn* 
calculated  in  Refs.  3  treat  the  fighter 
as  a  point  mass,  employ  a  and  bank  angle 
p  as  the  primary  “controls,”  and  employ 
steady  curves  for  the  coefficient*  C«_ 
and  Co.  Projected  on  a  vertical  plane. 
Fig.  1  reproduces  an  extreme  case  taken 
from  that  report.  The  vehicle  starts  and 
end*  at  the  same  point  in  space,  except 
that  it*  velocity  vector  &  longitudinal 
axis  are  exactly  reversed.  The  14JC  time 
advantage  thus  achieved  over  a  standard 
horizontal  turn*  with  the  same  initial 
and  final  conditions  but  constrained  a, 
can  probably  be  increased  by  the  clever 
use  of  som  unsteady  effects  examined  in 
Mhat  follow*.  Incidentally,  on*  of  the 
present  authors  published  a  rotational- 
dynamics  study  of  the  Ref.  3  maneuvers, 
in  which  it  was  shown  that  the  required 
aerodynamic  moment  control  is  feasible 
when  augmented  by  roughly  ±10*  thrust- 
vector  control  for  the  engine!*). 

In  a  timely  review  of  the  need*  and 
possibilities  relevant  to  the  design  of 
so-called  “agile”  aircraft,  Lang  and 
Francis4**  called  attention  to  likely 
significance  of  flow  unsteadiness  for 
enhancing  super maneuvers  (see  Figs.  6*  7 
and  8  of  Ref.  8).  The  present  paper 
undertakes  to  respond  to  their  call  by 
reexamining  cases  from  Refs.  3  in  the 
light  of  accurate  unsteady  «*asur«mint» 
that  have  just  begun  to  appear.  It  be¬ 
gin*  by  summarizing  som  of  the  recent 
f low-visual ization  tMts  and  aerod/naaic 
data  that  are  now  available.  Emphasis  is 
placed  on  those  tMts  wherein  transients 
were  performed  with  time  constant*  and 
rangM  of  a  which  CTosely  reproduce  the 
maneuvers  in  quastion.  The  most  suitable 
source  is  believed  to  be  a  prograe  con¬ 
ducted  by  one  of  the  present  authors  and 
partially  described  in  Refs.  9  and  10. 
Next  typical  exaeplM  from  Refs.  3  are 
recalculated  in  a  simplified  way  that 
permits  the  significance  of  unstMdy  ef¬ 
fects  to  be  assessed.  Finally,  and  in  a 
qualitative  effort  to  explain  the  root 
causM  of  t i me  lags  which  occur  in  the 
airloads,  results  of  an  empirical  "theo¬ 
ry"  are  coapared  with  lift  and  pitching 
moment  data  for  a  delta  wing  of  AR  *  1. 


Unsteady  Aerodynamics  ofL  low-AR  Hjpga 

With  regard  to  the  steady  flow  pat¬ 
terns  and  airloads  exDerienced  by  low-AR 
surfaces  or  complete  aircraft,  exposed 
to  low— speed  water  or  airflow  but  moder¬ 
ate  to  very  high  a,  the  literature  is 
extensive.  One  can  cite  general  surveys 
such  as  Refs.  11-13  and  the  forthcoming 
book  by  Rom41**.  Each  contains  useful 
theoretical  and  experimental  information 
with  an  emphasis  on  delta  planforms,  ei¬ 
ther  isolated  or  in  combination  with 
simple  bodies.  The  key  features  of  these 
flows  are,  of  course,  the  pattern  of  se¬ 
parated  vortices  which  exists  above  the 
lee  surfaces  and  the  manner  in  which  in¬ 
creasing  a  produces  progressive  develop¬ 
ment  of  instabilities.  Visualization  by 
the  illumination  of  smoke  traces,  etc., 
has  contributed  a  great  deal  to  their 
understanding.  A  seminal,  definitive  ex¬ 
ample  was  presented  by  Lowson4***  this 
year . 

Qualitative  work  on  the  consequen¬ 
ces  of  time— dependent  wing  motion  seems 
to  have  begun  in  Great  Britain  during 
the  early  1960's.  Thus  a  recent  investi¬ 
gation  by  Thompson,  Batill  &  Nelson4*** 
cites  Lowson’s  1964  discovery4 ***  of  a 
hysteresis  loop  in  the  location  of  L.-E. 
vortex  breakdown  or  "bursting”  above  a 
pitching  delta  model  with  sweep  n  -  80*. 
When  it  is  a  question  of  motions  which 
begin  to  reproduce  the  a-variation*  an¬ 
ticipated  for  super maneuvers,  however, 
sourcM  known  to  the  authors  are  limi¬ 
ted  in  both  numbers  and  scope.  Flow  vi¬ 
sualization  data,  focussed  on  the  beha¬ 
vior  of  vortex  breakdown  during  pitch 
transients  with  various  rangM  of  a,  are 
given  by  Nelson  and  coauthor*4 »**• 4 *■’ , 
Bad-el-Hak  6  Ho4  **» ■ ,mo> ,  Reynolds  &  Ab- 
tahi4***,  Atta  6  Rockwell ,mmt  and  Wolf- 
felt4--’ - 

Carrying  out  their  experiments  on  a 
delta  of  D-70*  and  a  wing-body  Soltani, 
Bragg  and  Brandon4***  present  force  and 
moment  coefficient*  for  siaple-harMnic 
pitching  between  a  »  0-  and  55°.  Their 
data  are  noteworthy  in  that  thrM  finite 
valuM  of  sldMlip  angle  $  are  included, 
as  are  unusually  high  valuM  of  Reynold* 
no.  Several  valuM  of  reduced  frequency 
k  —  over  the  range  of  really  practical 
InterMt  —  are  also  attained  with  the 
wing  alone.  The  other  Mjor  past  program 
involving  airload  Masurements  ha*  evi¬ 
dently  not  yet  received  full  publication 
but  certainly  deservM  citation.  In  the 
Netherlands  a  large  double-delta  Mdel 
with  ft’*  76*  Inboard  and  40*  outboard 
was  pitched  about  several  Man  m’s  and 
•t  aaplitudM  up  to  ±16*.  Som  data,  on 
time-dependent  surface  prMSurM,  normal 
force  and  pitching  Maent  are  given  by 
den  Boer  and  Cunningham4*-* ».  again  the 
Re’s  are  high,  ranging  froe  1.6  to  4.3 
million. 
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As  tool*  'for  examining  influences 
of  flow  unsteadiness  oh  high-ot  tnan*uv*r» 
th#  obvious  choice  must  b»  th*  airload 
measurements  reported  in  extenso  by  Jar- 
rah***'.  from  which  a  small  selection 
has  has  alrsady  appeared  in  Refs.  9  and 
10.  Six  components  of  fores  and  momsnt 
Here  taksn  by  strain-gaugs  balance  from 
sharp-edged  dslta  mods Is  of  AR’s  1.  l.S 
and  a  in  th*  SI  7  ft-by-10  ft  low-speed 
wind  tunnsl  at  NASA  Ames  Research  Ctr. 
Tn  thsss  tests  Aluminum  wings  mounted  on 
a  U-shaped  support  were  pitched  by  means 
of  hydraulic  actuation  about  an  axis  at 
two-thirds  midspan  chord.  Th*  angle  a 
was  varied  between  0”  and  values  up  to 
90" ,  either  in  a  ramp-like  fashion  or 
according  to  th*  sinusoid 

oc ( t )  *  <«mwx/2)C1  -  cosfltl  (1) 

Figure  2  is  a  schematic  of  the  apparatus 
u %ma  for  support  and  actuation.  Figure  3 
provide*  sketches  of  the  models,  inclu¬ 
ding  a  second  of  AR*1  which  was  used  for 
flow-visualization  tests  to  be  mentioned 
below.  The  reader  is  referred  to  Ref.  2 a 
for  extensive  details  on  orocedure,  data 
reduction,  error  estimation,  etc.  It  is 
added  that  only  four  of  the  six  airload 
components  can  be  reported,  because  side 
force  and  yawing  moment  were  always  zero 
within  the  accuracy  of  measurement.  The 
angle  0  was  held  to  zero  in  accordance 
with  th*  requirements  of  Refs.  1-3,  this 
being  the  constraint  used  to  prevent  de¬ 
parture  into  spins. 

Model  dimensions  are  given  in  mm  by 
Fig.  3.  Reynolds  nos.  based  on  midspan 
chord  e»  ranged  from  4.5  to  8.5  x  10*. 
The  dimensionless  parameter  characteri¬ 
zing  unsteadiness  for  both  sinusoids  and 
ramp  motions  is  chosen  to  be 

K  *  i««,*c./2V  <2) 

with  V  the  airspeed  and  the  time  deriva¬ 
tive  of  a  taken  at  Qt  »  w/2  <4ien  Eq.  <1) 
applies.  In  this  case  K  is  readily  con¬ 
verted  to  th*  more  conventional  reduced 
frequency  kj  multiplication  would  be  by 
2/i»  when  a***  -  90*.  Values  of  K  fro*  0 
(steady  flow)  up  to  0.08  were  obtained 
in  th*  wind  tunnel.  This  may  be  compared 
with  s  msximu*  of  about  K  *  0.1  for  the 
maneuver*  analyzed  in  Ref*.  3,  where  in 
most  case*  the  parameter  was  less  than 
0.05.  — 

Th#  new  data  are  shown  in  Ref -  26 
to  correlate  quit*  satisfactorily  with 
steady-flow  counterparts  over  th#  range* 
of  a,  R*  and  Afi  where  th#  latter  prove 
to  be  available.  Comparison*  are  also 
made*  where  possible,  with  airloads  from 
Ref*.  23.  24  and  25  »  again,  systematic 
and  unexplained  discrepancies  are  not 
found.  Figures  566  typify  the  present 
steady-flow  aerodynamic  coefficients.  A* 
in  other  examples  which  follow,  plots  v* 
a  are  shown  as  continuous  curve*  since 


th*  data-reduction  system  stored  infor¬ 
mation  at  interval*  of  less  than  one  de¬ 
gree.  One  exception  is  Fig.  6,  wherein 
lift  measurements  appear  as  squares  for 
the  AR»1.5  model  and  are  compared  with 
two  other  experimental  tource*1*'’1  •  <*•> 
and  with  the  theory  of  Polhamum‘"*» .  The 
wing  of  AR“2  is  selected  for  most  of  th* 
illustrations  her*  because  this  is  the 
on*  later  used  in  maneuver  analyses.  All 
plotted  coefficients  are  defined  accoi — 
ding  to  standard  aeronautical  practice. 
For  example,  lift,  pitching  moment  and 
rolling  moment  are,  respectively, 

Cu  >  L/(<r/2>V«^  C„  »  M/<oV2)V*c^ 

Cn  -  R/(<r/2)V*bS’  <3a,b,c) 

Moment  R  acts  about  the  midspan  axis. 
Pitching  moment  fi  is  positive  nose-up 
about  th*  77X-chord  axis,  so  as  to  en¬ 
sure  that  th*  values  remain  generally 
positive,  but  th*  reader  is  reminded 
that  the  pitch  axis  is  at  two-thirds  the 
midspan  chord  for  unsteady  testing. 

For  th*  AR»2  delta,  Figs.  6  throuqh 
9  show  the  histories  of  five  aerodynamic 
coefficients  as  oi  is  varied  sinusoidally 
through  on*  cycle  from  0"  to  90"  and  re¬ 
turn.  Arrows  on  th*  curve*  give  the  di¬ 
rection  of  motion.  Th#  pitch-rat*  para¬ 
meter  K  —  taken  as  the  measure  of  flow 
unsteadiness  —  increase*  from  0.01  up 
to  th*  intermediate  value  0.04  when  on* 
goes  through  these  four  figures. 

.  Fro*  careful  study,  a  number  of  ob¬ 
vious  conclusions  can  be  drawn,  nearly 
all  of  which  apply  for  three  models  and 
for  both  th*  sinusoidal  and  raep  test* 
that  went  wall  past  maximum  lift.  Even 
at  values  of  K  below  0.01  lift,  normal 
force  and  drag  significantly  exceed  th* 
corresponding  steady  value*  at  a’s  above 
20-25"  when  this  angle  is  increasing  but 
fall  well  below  steedy-flow  on  th*  down- 
stroke.  This  overshoot  becomes  larger  as 
K  increase*.  Its  remarkable  magnitude  is 
estiaated,  for  example,  by  comparing  th* 
curve*  of  normal  force  between  Figs.  9  &■ 
4.  Th#  peak  of  th*  graph  move*  to  higher 
a  and  at  K-0.04  exceeds  it*  steady-state 
value  by  over  SOX .  A*  can  be  concluded 
fro*  prior  tests  and  free  th#  flow-visu¬ 
alization  analyses  in  Ref.  26,  behavior 
of  thia  sort  is  connected  with  delay*  in 
th#  breakdown  or  "bursting"  of  th#  L.-E. 
vortex  system  on  th*  upstroke,  followed 
by  a  lag  in  its  reestablishment  as  th*  a 
returns  toward  zero. 

Th#  time  histories  of  pitching  mo¬ 
ment  reveal  the  same  increasing  trends, 
and  it  can  be  inferred  that  th*  center 
of  normal  force  move*  forward  on  th* 
wing  as  th*  chordwis*  location  of  vortex 
breakdown  proceeds  forward  during  th* 
upstroke.  It  is  regarded  a*  important 
for  th#  feasibility  of  high-«  maneuver* 
that,  unlike  what  has  been  observed  on 
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pointed  bodies  of  revolution  end  delta 
Mings  in  sideslip*  rolling  moments  stay 
consistently  very  small.  At  no  time  is 
the  center  of  lift  found  to  move  off  the 
Ming  centerline  by  more  than  about  0.5% 
of  the  wingspan.  These  same  observations 
hold  for  the  AR“1 .5  model.  At  AR»1  roll¬ 
ing  moments  show  a  more  erratic  behavior 
between  a  ■  23*  and  5^ »  especially  at 
very  low  K.  One  believes*  however*  that 
their  excursions  are  not  beyond  the  abi¬ 
lity  of  aerodynamic  controls  to  balance. 

Three  additional  figures  are  inclu¬ 
ded  as  representative  of  the  extensive 
data  collected  during  this  program.  In 
Fig.  10  are  plotted  the  five  aerodynamic 
coefficients  for  the  AR» 1 . 3  model  at  the 
high  K«0.06.  The  mild  oscillations  seen* 
for  instance*  in  all  the  curves  during 
downstroke  are  not  regarded  as  indica¬ 
tions  of  experimental  inaccuracy.  They 
are  reproducible  in  repeated  tests  and 
are*  therefore*  in  need  of  explanation. 
Figure  11  demonstrates  the  influence  of 
parameter  K  on  normal  force  for  AR“1. 
The  solid  curves  here  are  for  very  slow 
variation  of  a|  it  is  not  known  whether 
the  small  differences  between  up-  and 
downstroke  constitute  soee  sort  of  hys¬ 
teresis  or  merely  test  imprecision.  The 
final  example*  Fig.  12*  shows  the  effect 
on  normal  force  of  varying  Re  between 
.45  and  .85  million  for  the  AR«1.3  wing 
at  100.02.  As  in  other  experiments  that 
have  been  conducted  on  deltas  with  sharp 
leading  edges  over  considerably  wider 
ranges  of  this  parameter*  it  is  not  felt 
that  any  significant  influence  or  Re  on 
resultant  airloads  can  be  detected. 


Unsteady  Effects  on  Turning 

The  measurement  program  reviewed  in 
the  preceding  section  provides*  perhaps 
for  the  first  time*  a  chance  to  quantify 
the  potential  for  enhanced  fighter  agi¬ 
lity  inherent  in  the  remarkable  flow  un¬ 
steadiness  over  pitching  delta  wings. 
For  many  years  the  favorable  and  unfa¬ 
vorable  effects  of  "dynamic  stall"  have- 
been  studied  for  wings  of  moderate  to' 
high  AR>  rotors*  wind  turbines*  etc.  It 
is  foreseen  that  very  detailed  analyses 
of  this  subject  for  low— AR  aircraft*  in¬ 
cluding  trajectory  optimisations  and 
combat  simulations*  will  be  required  be¬ 
fore  new  designs  and  operational  proce¬ 
dures  can  be  adopted.  Certainly  unsteady 
wind-tunnel  testing  of  complete  models 
will  become  routine  practice.  At  the 
current  level  of  understanding*  however* 
a  much  simpler  approach  seems  all  that 
is  justified. 

For  the  present  investigation*  it 
was  therefore  decided  simply  to  reana¬ 
lyse  the  response  of  the  "generic"  all — 
craft  of  Ref.  7*  as  it  executes  turning 
maneuvers  defined  according  to  the  time 
histories  of  the  controls  a(t>  and  #(t> 
taken  from  Well  et  a 1 . ’ ■ * .  Based  on  the 


properties  of  typical  fighters  dis¬ 
cussed  in  Refs.  3  and  R»n*o«'*°,i  this 
vehicle  has  a  mass  of  10*617  kg  and  coi — 
responding  moments  S-  product  of  inertia. 
The  double-delta  wing  of  Ref.  7,  with  S« 
37.7  m>*  is  replaced  by  a  single  sharp- 
edged  delta  of  AR-2  so  as  to  permit  di¬ 
rect  use  of  data  like  that  on  Figs.  6-9. 

As  in  Refs.  3*  it  is  assumed  that 
rotational  motions  can  be  accomplished 
so  rapidly  as  to  place  no  constraint  on 
the  vehicle  dynamics.  Accordingly*  tra¬ 
jectory  computations  can  be  carried  out 
as  if  it  were  a  point  mass  m.  Maneuvers 
are  then  described  by  the  time  histories 
of  airspeed  V*  flight-path  angle  r  above 
the  horizontal*  and  azimuth  angle  X  — 
measured  from  a  horizontal  datum  clock¬ 
wise  around  to  the  vertical  plane  which 
contains  V.  The  state  vector  V*  r,  X  is 
governed  by  three  nonlinear*  first-order 
differential  equations*  as  followst 

mV  -  Tcosa  -  D  -  mgsinT  (4) 

mVT  -  CTsina  +  Llcosd  -  mgcosr  (5) 
(mVeosDX  «  CTsina  *  Lisin**  (6) 

There  are  three  auxiliary  kinematic 
relations  for  rate  of  change  of  altitude 


h  and  two  horizontal  coordinates  xs* 
in  a  (no— wind)  earth— fixed  triad* 

y» 

• 

*• 

«  V  cosT  cosX 

(7) 

V- 

•  V  cosT  sinX 

(8) 

* 

-  V  sinT 

(9) 

Equations  (7>-<9>  can  readily  be  used  to 
construct  the  trajectory  in  space*  but 
results  of  this  sort  at  not  given  here. 

Before  presenting  some  solutions  of 
the  system  (4>-C6>,  a  few  remarks  are  in 
order.  In  Refs.  3  a  fourth  equation  was 
discussed  which  connects  the  rate  of  de¬ 
crease  of  mass  m  to  thrust  T  and  engine 
fuel-consumption  data.  All  the  maneuvers 
of  interest  here  occur  in  such  short  in¬ 
tervals*  however*  that  m  is  essentially 
constant.  There  is  a  body  axis  x  along 
the  zero-lift  direction  and  inclined  at 
angle  a  above  the  flight  path*  but  it  is 
bank  angle  about  velocity  V.  positive  to 
depress  the  right  wing  below  .the 
horizontal.  In  Refs.  3*  T  and  a  drag- 
brake  rotation  angle  are’ used  as  auxi¬ 
liary  controls.  Along  with  a  and  d*  the 
values  used  for  these  are  taken  straight 
froe  that  source. 

Except  for  the  coefficients  C*_  and 
Co*  all  inforeation  needed  for  the  tra¬ 
jectory  calculations  can  be  taken  from 
the  large  appendix  of  the  DFVLR  report* 
part  of  Refs.  3  and  supplied  to  the  au¬ 
thors  by  Dr.  Well.  As  a  quantitative  ap¬ 
proach  to  the  primary  objective  of  this 
paper*  a  scheme  has  been  devised  to  make 
direct  comparisons  between  similar  tra- 
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4  jtctorie*  determined*  respectively*  from 

quasi-steady  aerodynamic  information  and 
corresponding  unsteady  data.  Steady  and 
unsteady  coefficients  were  taken  fom 
the  same  value  of  a  at  each  value  of  t 
called  out  in  a  numerical  integration  or 
Eqs.  (4) -(6).  The  steady  C*_  and  Co  come 
from  the  dotted  and  dash-dot  curves  on 
^lg.  4*  respectively.  The  unsteady  data 
come  from  curves  like  :ho3B  or.  Figs.  6-9 
with  the  value  of  K  estimated  as  closely 
as  pcssible  from  the  prescribed  history 
of  the  control  <*.  Obviously*  the  upper 
branch  cf  the  curve  is  used  when  *  is 
increasing  and  the  lower  when  ut  is  de- 
creasi ng . 

Since  the-  lift  and  rirag  information 
employed  here  does  not  agree  exactly 
with  that  i.i  Rvfs.  3  Ccf.  Fig.  1  of  the 
paper  in  J.  Guidance.  Control  and  Dvn.)* 
one  muse  ensure  that  the  tra jector ies 
calculated  from  steady  data  are  reison- 
ab i v  close  tc  Die  arjther.  Comparisons 
are  made  in  certain  of  the  following  ex¬ 
amples  *  The  approach  used  here  in  p-der 
to  isc  la's  unsteady  effects  is  believzd, 
however,  to  t-e  the  only  logical  one. 

u’-e  maneuvers  chosen  for  study  ..re 
dimple*  a  very  recent  article**’ *  shows 
.hat  they  rio  resemble  several  of  those 
being  used  tor  Fi  ight  demonstration  in  a 
rrorjrain  cohd-icted  by  USAF  and  NASA.  Most 
cases,  as  in  earlier  analytical  studies, 
emphasize  a  reor  )  ente.t  ton  uf  the  fuse¬ 
lage  axis  and/or  the  velocity  /-acto.  in 
minimum  time  from  a  given  i  dial  stale. 
It  is  generally  agreed  that  these  objec¬ 
tives  are  closely  associated  with  aa>,  i- 
mi m  attainable  values  of  the  pitching 
angular  velocity  D.  This  quantify  is* 
therefore.  the  figure  of  merit  employed 
here.  No  .attempt  ; «  made  to  meet  pre¬ 
scribed  final  conditions  or  to  optimize* 
since  such  sophistication  is  beyond  whr.t 
can  be  justified  in  the  light  of  present 
approximations.  Results  of  the  selected 
examples  are  new  listed  and  illustrated* 
each  case  fceinq  identified  with  its  num¬ 
ber  from  Refs.  3. 

(1)  4. 2. 2-1  --  A  horizontal  turn,  with 
the  objective  of  rotating  the  velocity 
vector  through  180“.  Initial  velocity  is 
100  m/s,  and  <x  is  constrained  to  be  less 
than  about  28*.  This  is  clearly  not  a 
" supermaneuver " -but  is  used  to  provide  a 
standard  of  reference.  Figure  13  plots 
the  prescribed  angles  of  attack  and  bank 
as  spline  curve  fits  to  data  giver  for 
six  time  instants  on  page  6  of  the  Refs. 
3  appendix.  Fiqures  14*  IS  and  16  show* 
respectively*  the  histories  of  airspeed* 
angular  velocity  and  normal  load  factor 
n«  over  the  nearly  8  s  required  for  the 
conventional  turn.  Note  that  the  curves 
computed  with  steady  (dotted)  &  unsteady 
(solid)  airloads  agree  closely*  as  anti¬ 
cipated  in  view  of  the  closeness  of  the 
coefficient  plots  (e.g.*  Fig.  11)  below 
a»88".  (The  scales  are  greatly  expanded 
on  Figs.  14-16.)  Incidentally,  n«  is 
calculated  in  conventional  fashion,  e.g. 


from  Eqs.  (12)  of  Refs.  3.  Convergence 
studies  using  progressively  lower  steps 
in  the  time  integrations  have  shewn  all 
results  to  be  accurate  within  the  pre¬ 
cision  of  plotting.  Similar  good  agree¬ 
ment  was  obtained  here  with  the  varia¬ 
tions  of  n„,  heading  and  flight-path 
angles.  It  can  be  concluded  that  flow 
unsteadiness  does  not  improve  the  exe¬ 
cution  of  maneuver  4. 2.9-1. 

(2)  4. 2. 6-1  —  A  vertically  upward  turn 
whose  objective  is  to  reverse  the  axis 
of  the  fuselage  but  with  no  final  con¬ 
straint  on  oi  or  flight-path  angle.  The 
initial  velocit /  is  100  m/s*  again  se¬ 
lected  because  high-a  agility  is  best 
demonstrated  at  low  speeds*  where  ex¬ 
cessive  acceleration  can  be  avoided. 
The  "controls"  -ire  plotted  in  F;q.  17 
over  the  5.4  s  required  for  the  mareu- 
/er  according  to  Refs.  3.  The  wings  re¬ 
main  essentially  level*  and  it  is  re¬ 
marked  that  this  variation  of  a  seems 
more  realistic  than  the  fluctuating  one 
shown  on  page  lo2  of  thB  Refs.  3  appen¬ 
dix.  From  Fig.  17  and  the  airspeed  his¬ 
tory  on  Fig.  18*  it  was  determined  that 
the  best  estimates  of  unsteady  airloads 
should  be  taken  from  upward-ramp  tests 
at  K  *  0.03.  Given  these  data,  Eqs.  (*♦/  — 
(6)  yield  the  responses  graphed  on  Figs. 
13-21.  The  increased  drag  at  high  a  pro¬ 
duces  the  small  unsteady  reduction  in  0 
on  ciq.  18.  Pitch  rates  on  Fig.  19  are 
seen  to  be  considerably  higher  with  un¬ 
steadiness  accounled  f c  ■  *  especially  In¬ 
ward  the  end  of  the  turn.  Since  the  path 
angle  T  (Fig.  20>  is  one  possible  dete;  - 
minant  of  when  success  is  achieved.  one 
sees  that  a  value  of  9>»  is  reached  0.4 
s  (or  almost  10%)  Faster  than  prior  es¬ 
timates  would  indicate.  Figures  IE  &.  20 
contain  triangular  points  taken  right 
from  the  tables  of  Refs.  3  and  suggest* 
m  this  case*  that  the  differences  bet¬ 
ween  the  airloads  used  thei  e  and  here  do 
not  lead  to  large  discrepancies.  Figure 
21*  finally*  implies  that  the  normal  <c- 
celerctions  and  the  associated  struct¬ 
ural  loads  are  not  affected  unfavorably 
by  unsteadiness  and*  in  Fact,  are  some- 
wnat  lower  near  the  end  of  the  turn. 

(3)  4. 2. 7-1  --  A  ve'-t  ica  1  J  /  upward  turn 
whose  objective  is  to  reverse  directions 
of  both  the  fuselage  axis  and  airspeed 
vector  at  the  top.  Again.  initial  V  is 
100  m/s  and  a  is  unconstrained.  It  is 
worth  mentioning  that  no  time  allowance 
is  made  at  the  end  of  these  two  maneu¬ 
vers  for  the  180“-roll  required  for 
bringing  the  cockpit  upright.  as  in  an 
Immelmann  turn.  The  assumption  is  that 
these  adjustments  can  be  made  rapidlvi 
thev  are  unlikely  to  benefit  much  f,  om 
unsteady  flow.  The  controls  for  4. 2. 7-1 
appear  on  Fig  22.  Figuras  23-24  graph 
corresponding  histories  of  V,  Q  and  n, . 
In  this  case  the  11  -  cosfltj  unsteady 
data  for  K  ■  0.04  were  found  to  provide 
the  best  approximations.  As  with  4. 2. 6-1 
unsteady  effects  reduce  the  airspeed  and 
yield  (here  mere  modest)  improvements  .n 
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turn  rat*  during  the  high-a  portion  of 
the  turn.  Load  factors  (Fig.  35)  prove 
so«snh*t  higher  up  to  about  t  »  4  s ,  but 
this  is  not  believed  to  be  an  unsteady 
effect  because  the  a’s  are  relatively 
low  in  this  range. 

(4)  4.3—2  —  As  sketched  in  Fig.  1,  this 
is  a  rather  violent  supermaneuver  where¬ 
in  the  aircraft  starts  and  ends  at  the 
same  point  in  space  but  with  both  the 
axis  and  airspeed  vector  reversed.  Ini¬ 
tial  speed  is  100  m/s,  and  a  is  uncon¬ 
strained.  In  Refs.  3  the  estimate  is 
that  a  14X  time  advantage  results  from 
turning  this  way  rather  than  using  banks 
in  a  horizontal  plane  under  constrained 
at*  even  with  both  maneuvers  performed 
optimally.  It  is  evident  from  Fig.  1* 
however ,  that  automatic  stabilization 
and  control  will  be  required  to  hold 
zero  sideslip  and  otherwise  follow  this 
sinuous  path. 

Figure  36  furnishes  curve-fits  to 
the  at ( t )  and  P(t>  from  Refs.  3.  Again 
the  sinusoidal  airload  data  for  K  -  0.04 
seemed  most  suitable  for  supplying  un¬ 
steady  effects  during  the  high-or  transi¬ 
ent.  Figures  37-39  give  the  steady-vs.- 
unsteady  comparisons.  On  the  airspeed 
curves*  Fig.  37,  the  triangles  show  that 
present  steady  results  almost  coincide 
with  those  of  Refs.  3*  except  quite  near 
to  the  maneuver’s  end.  Figure  28  yields 
the  interesting  information  that  there 
are  substantial  unsteady  improvements  in 
turn  rat*  through  the  increasing-*  phase 
but  that  due  to  combinations  of  dynamic 
•*f*cts  nothing  is  lost  during  the  down- 
stroke.  Load  factors  (Fig.  29)  show  some 
modest  increase  during  the  first  sharp 
pull-up  but  remain  slightly  below  their 
steady  values  for  the  rest  of  the  turn. 

(3)  4. 2. 7-8  —  A  reversing  vertically 
upward  turn  similar  to  4. 2. 7-1,  but  with 
initial  velocity  300  m/s  and  a  limited 
below  70- .  Several  other  of  the  Refs.  3 
maneuvers  have  been  analyzed  in  a  manner 
similar  to  the  above,  and  it  would  be 
misleading  to  imply  that,  in  all  cases, 
the  influence  of  unsteadiness  will  give 
rise  to  faster  turns  with  no  penalty  in 
terms  of  structural  loads.  In  general, 
it  is  found  that  examples  starting  above 
V  -  lOO  m/s  are  not  so  favorable.  The 
low-altitude  "corner  velocity”  for  the 
aircraft  studied  here  is  around  143  m/s, 
and  on*  speculates  JEftat  beyond  this  air¬ 
speed  there  is  little  to  be  gained,  be¬ 
cause  of  load-factor  limitations. 

Figures  30-33  plot  the  same  infor¬ 
mation  for  4.E.7-8  as  has  been  given  for 
the  foregoing  cases.  The  sinusoidal  air¬ 
load  data  for  K  -  0.08  were  used  in  the 
unsteady  calculations.  Figure  31  shows 
that  drag  overshoot  causes  some  bleeding 
off  of  airspeed,  but  the  Fig.  38  turn- 
rate  advantages,  which  occur  only  below 
the  "corner,"  are  not  impressive.  The  n« 


peak  on  Fig.  33  is  unacceptable*  and  un¬ 
steadiness  is  seen  to  increase  it  by  ap¬ 
proximately  one  “g." 


A  Qualitative  Aerodynamic  Theory 

There  is  a  long  history  of  attempts 
to  extend  to  time-dependent  motions  the 
many  steady-flow  theories  that  have  been 
proposed  to  account  for  the  organized 
pattern  of  free  vortices  that  develops 
above  slender,  pointed  wings  as  a  is  in¬ 
creased  (Roa*1*'  summarizes  the  latter 
thoroughly).  The  contributions  of  Low— 
son**-*,  Dore****  and  Randall****’  from 
England  in  the  196C>’s  deserve  first  men¬ 
tion.  More  recently  adaptations  of  panel 
methods  have  been  published  (e.g..  Levin 
&  Katz****),  and  van  Niekerk***’  modi¬ 
fied  the  leading-edge-singularity  scheme 
of  Polhamus***’  to  account  I’or  time- 
varying  a  on  a  delta  wing.  Two  comments 
are  offered  regarding  all  ,f  this  work. 
The  first  is  that,  for  slues  of  para¬ 
meters  K  or  k  typical  of  supermaneuvers • 
the  airloads  they  predict  are  essenti¬ 
ally  quasi -steady .  Secondl*-,  none  seems 
capable  of  modelling  vortex  breakdown. 
Indeed,  the  authors  believe  that  a  fully 
rational  theory  would  have  to  be  based 
on  the  methods  of  computational  fluid 
dynamics  and  would  have  to  account  for 
large  volumes  of  separated,  turbulent 
flow.  Analyses  of  this  sort,  feasible  of 
extension  to  a’s  as  large  as  90-,  are 
probably  well  beyond  the  scope  of  even 
the  most  powerful  current  CFD  methods. 

Granted  the  impossibility  of  relia¬ 
ble  predictive  tools  in  the  near  future, 
on*  is  forced  to  conclude  that  unsteady 
wind-tunnel  testing  is  the  only  alterna¬ 
tive  available  to  the  designers  of  agile 
aircraft.  Every  flight  vehicle  must  pass 
through  a  preliminary  design  phase,  how¬ 
ever,  when  its  configuration  is  not  well 
enough  established  to  permit  model  con¬ 
struction.  On*  is  perhaps  justified,  in 
such  a  situation,  when  he  puts  forth  a 
purely  "empirical”  or  "qualitative"  at¬ 
tempt  to  reproduce  the  principal  fea¬ 
tures  of  a  phenomenon. 

Any  such  approach  must  rely,  first 
of  all,  on  estimates  of  vortex  breakdown 
whose  hysteresis  is  known*1*’ •*•» »  to  be 
the  controlling  cause  of  unsteadiness. 
At  the  higher  Reynolds  nos.  breakdown  is 
a  rather  sudden  process.  Let  is  be  as¬ 
sumed  that  information  is  available  froa 
sources  like  Refs  ,  16,  18—84  and  36  on  a 
quantity  »•,  as  it  varies  with  *  during 
prescribed  motions  (ramp,  sinusoid)  for 
a  useful  rang*  of  a  parameter  like  K.  x* 
is  her*  defined  to  be  the  forward  limit 
of  an  identifiable  breakdown  region. 

As  an  examol*  of  the  sort  of  data 
needed.  Fig.  34  from  Jari*ah«**»  gives 
estimates  of  the  angles  where  x*  passe* 

the  73X— midspan— chord  station  of  an  AR-i 
delta  on  the  up-  and  down-strokes  of 
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sinusoids  between  a  -  0-  and  90.  From 
flow  visualizations  this  quantity  proved 
to  depend  on  K  only?  for  the  upstroke, 
it  settles  down  fairly  quickly  to  values 
near  45*.  Jarrah****  provides  similar 
data  for  chordwise  stations  at  50V.  and 
the  trailing  edge.  The  results  are  con¬ 
sistent  between  sin.  aoidal  and  ramp  mo¬ 
tions,  and  they  are  found  to  agree  with¬ 
in  measurement  accuracy  with  the  data  of 
Reynolds  and  Abtahi<»»>  for  an  almost 
identical  wing  model. 

In  the  spirit  of  Ericsson’s  imagi¬ 
native  insights  (Ref.  37  contains  recent 
examples),  let  it  be  assumed  that  the 
distribution  of  aerodynamic  force  per 
unit  chordwise  distance  is  made  up  of 
three  parts:  <1)  A  portion  determined 
from  the  rate  of  change  of  crossflow 
momentum,  in  a  manner  resembling  the 
low-AR  theory  of  Jones *»•»  but  with  the 
slabs  of  fluid  taken  normal  to  the  wing 
surface  at  a.  (2)  A  portion  calculated 
on  a  quasi-steady  basis  by  rotating  the 
L.-E.  suction  force  through  90“ ,  as 
proposed  by  by  Polham.is'*” .  At  higher 
<*’*»  loads  (1)  8,  (t)  act  only  ahead  of 
xwt<x(t)3.  (3).  A  portion,  dominant  at 
the  higher  a’s  and  calculated  from  the 
Betz**1**  crossflow-drag  model.  This 
turnr  out  nearly  proportional  to 
Co*sin*a.  C  itself  is  chosen 

empirically  from  the  drag  at  <*»90“ 
measured  for  a  given  K . 

For  pitching  «(t>  about  a  fixed  ax¬ 
is  at  two-thirds  midspan  chord,  formulas 
for  normal  force  and  moment  about  that 
same  axis  are  given  in  Eqs.  (10)  &  (11), 
which  follow.  The  moment  can  be  trans¬ 
ferred  to  77X— chord  in  the  usual  way  for 
experimental  comparisons.  It  seems  con¬ 
sistent  to  assume  that  the  force  resul¬ 
tant  acts  perpendicular  to  the  wing,  so 
that  lift  and  drag  are  just  the  cosa  and 
sina  components  of  C„.  In  these  formulas 
the  three  portions  listed  above  are  con¬ 
tained  sequentially  in  the  three  braces, 
"sin"  &  "cos"  are  abbreviated  “s"  6  "c." 

Cm  -  <»/6>  (  x»/c.)  *ctnrK6s2o-B(4c./V>c<* 
+16(4c./V)  ( x»/c.  )ca+(9c.*/V«  >C3(x»/c.  >  * 
-(8x»/3c.)J>  «■  («/6s0)  (xa/c.)  *<6s*a 
-►8 <  4c.  / V )  C  <  x./c.  ) -1  Isa 
+  (4c,./V)*C3<x./c.>*-#-(B/3)  -  (16x«/3c.)J> 
+  Ca,){»««+(l/lB)(ic./V)>)  (10) 

Cm.  -  -8e<  Xw/e.  >  *  <e./C>ctnTKC  < 

-1 3se<rM*c./V>  C3<x.Zc  >  *-4<x./c.  >+4/33c« 
♦<lfe.«/V*)C<3/5)  <xw/e.)«-<x»/c.)* 
♦(4x,/9c.)3>  -  ( w/6sT)> ( Xm/c.) I ( c./E )* 
*<4C<x,/c.)-Us«a+  <4c./VK3<Xe/c.)« 

-(  16x»/3c.)+(B/3>  Isa-Mdc./V)  *C-<  16/9) 

♦ (lAx./Sc. ) -6 ( x./c. ) *♦ < I e/S > < x ,/c . ) > 

-Cm.  (Cm/eX  (1/9)  < 4c ./Visa 

-(  1/135)  (4c./V>  *  >  (ID 

For  the  delta  wing  of  AR-1 ,  Figs. 
33  and  36  show  as  solid  curves  the  pre¬ 
dictions  of  Eqs.  ( 10)  &  (11),  plotted 
vs.  «  for  K-0.06.  The  moments  of  Fig.  36 
are  about  the  77X -chord  axis,  so  that  •. 
direct  comparison  is  possible  with  the 


measured  data,  plotted  as  dash-dots. 
This  aspect  ratio  was  selected  because 
spline  fits  could  be  used  for  Xm  esti¬ 
mates  drawn  from  Reynolds  &  Abtahi*"1* 
and  the  Ref.  26  flow  visualizations.  All 
that  can  be  stated  is  that  this  first 
attempt  at  a  "theory"  reproduces  the 
qualitative  behavior  Df  the  airloads. 
Similar  calculations  for  AR’s  1.5  and  2 
yisld  comparable  or  better  accuracy.  It 
also  appears  possible,  by  more  realistic 
handling  of  the  parameter  x„  to  get 
better  agreement  over  the  high-oc  por¬ 
tions  of  these  curves.  Finally,  resort 
to  more  precise  theories  like  that  of 
Dore ‘ ■■ *  is  likely  to  give  improve— 
m»r>ts  in  the  a-range  where  breakdown  has 
little  influence  on  the  loading. 

Concluding  Remarks 

A  program  has  been  reviewed  of  air¬ 
load  measurements  on  a  family  of  low-Arl 
delta  wings  with  sharp  leading  edges, 
subjected  in  the  wind  tunnel  to  large- 
amplitude  pitch  transients  involving  g’s 
going  as  high  as  90.  Rather  small  val¬ 
ues  of  the  pitch-rate  parameter  K  were 
u»md«  representative  of  maneuvers  anti¬ 
cipated  for  "agile”  aircraft.  Even  for 
»se  modestly  unsteady  motions,  it  is 
found  that  force  and  moment  overshoots 
can  exceed  by  50%  their  steady-state 
counterparts.  The  explanation  lies  in 
the  hysteretic  behavior  of  the  break¬ 
down  location  of  L.-E.  vortices. 

By  means  of  examples  based  on  low- 
speed,  high-a  maneuvers  from  the  litera¬ 
ture,  an  attempt  is  made  to  demonstrate 
that  considerably  higher  turn  rates  can 
be  achieved  than  would  be  predicted  from 
steady-state  airload  data.  This  enhanced 
maneuver abi 1 i ty  is,  by  no  means,  always 
accompanied  by  a  penalty  in  tarns  of 
load  factor  and  associated  pilot  discom¬ 
fort  or  structural  ovarstresa.  Not  all 
tha  cases  studied  show  these  advantages, 
however,  because  the  ieproved  agility 
appears  to  exist  only  at  flight  speeds 
well  below  the  "corner  velocity."  It  is 
believed  that  designers  of  these  vehi¬ 
cles  (and  their  control  systems)  should 
certainly  saek  to  take  account  of  the 
potentialities  of  flow  unsteadiness. 

The  paper  concludes  by  proposing  a 
very  approximate  theoretical  model  which 
tries  to  include  the  breakdown  hystere¬ 
sis  as  part  of  a  three-term  representa¬ 
tion  of  the  unsteady  chordwise  load  dis¬ 
tribution.  The  resulting  estimates  for 
normal  force  and  moment  due  to  pitching 
eotion  exhibit  the  same  features  found 
in  test  data,  but  more  refinement  will 
obviously  be  needed  before  this  model 
has  any  chance  of  quantitative  success. 
It  is  put  forth  in  the  conviction  that  a 
wholly  rational  theory  eust  await  ex¬ 
tensive  developments  in  the  field  of 
CFO.  Lacking  such  tools,  however,  one 
concludes  that  wind-tunnel  tests  on  pit¬ 
ching  models  of  agile-aircraft  designs 
will  furnish  the  only  pre-flight  source 
of  the  information  required  to  analyse 
♦■heir  maneuvers. 
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Fig.  1.  Vehicle  orientations  and  trajec¬ 
tory  in  vertical  plane  for  supermaneuver 
.involving  reversal  of  velocity  and  fuse¬ 
lage  axis.  Initial  and  final  positions 
are  at  the  same  point  in  space. 


in  the  wind  tunnel . 
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Fig.  3.  Models  used  for  airload  tests 
and  flow  visualization. 
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4.  Plot  of  four  coefficients  vs.  c 
steady  flow  oast  delta  model  of 
Re/nr  Ids  no.  Re  *  5-9x10*. 
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~ig.  5.  Steady  lift  coefficient'  of  f>R«c 
delta  model.  compared  with  theory  of 
Poihamus'*”  and  tests  of  Lawford  •m'r*  S, 
Bert i et  t 1 ** *  . 
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Fig.  6.  Coefficients  of  normal  force, 
drag,  rolling  moment  (essentially  zero), 
{lift  and  pitching  moment  about  77%-chord 
;a;:is>  these  are  plotted  vs.  or  for  AR«3 
delta  performing  the  Cl  -  cosOtl  pitch 
.maneuver.  Unsteadiness  parameter  is  K« 
0.01.  and  Re  -  4.5x10*. 
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.  7.  Same  as  Fig.  6  for  K»0. 02  and 
.4x10*. 
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.  S.  Same  as  Fig.  6  for  K-O.03  and 
.5x10*. 
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g.  9.  Sam*  as  Fig.  6  for  K-0.04  and  R« 
4. 5x10*. 


Fig. IS.  Effect  of  Reynolds  no.  on  normal 
fores  history  for  AR-1.5  delta  it  K  > 
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Fig.  10.  Same  as  Fig.  6  for  AR“1.S  delta 
at  the  high  value  K— O . OA . 
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Fig.  13.  Time  histories  of  a  and  ft  a* 
“controls"  for  minimum-time  conventional 
turn  (No.  4.2.2— 1  of  Ref*.  3). 
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Fig.  11.  Normal  force  vs.  a  for  AR- 1 
elta  performing  the  Cl  -  cosQtl  pitch 
tneuver  at  four  values  of  K. 


Fig.  14.  Histories  of  airspeed  V,  calcu¬ 
lated  for  Maneuver  4. 2. 2-1  with  steady 
(dashed)  and  unsteady  (solid  line)  air¬ 
loads. 
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Fig.  80.  Histories  of  flight-path  angle 
r  for  Maneuver  4. 8. 6—1.  Triangles  .are 
point*  found  with  data  fro*  Ref*.  3. 
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Fig.  21.  Historian  of  normal  load  factor 
n_  for  Maneuver  4. 2. 6-1. 
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Fig.  24.  Historian  of  pitch  rate  C  for 
Maneuver  4.2.7- 1. 
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Fig.  22.  Time  histories  of  a  and  d  as 
"controls"  for  minimum-time  vertically- 
upwar d  turn  without  a— constraint  but 
with  reversal  of  velocity  and  fuselage 
axis  (No.  4,2. 7-1  of  Refs.  3). 
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Fig.  23.  Histories  of  normal  load  factor 
n.  for  Maneuver  4. 2. 7-1. 
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Fig.  23.  Histories  of  airspeed  V»  calcu¬ 
lated  for  Maneuver  4. 2. 7-1  with  steady 
(dashed)  and  unsteady  (solid  line)  aii — 
loads. 
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Fig.  26.  Time  histones  of  a  and  d  an 
“controls"  for  minimum-time  reversal 
without  a-constraint  (Fig.  1  and  No. 
4.3-2  of  Refs.  3). 
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